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The ability to simulate one Hamiltonian with another is 
an important primitive in quantum information processing. 
In this paper, a simulation method for a rbitrary a z (g) a z in - 
teraction based on Hadamard matrices ( quant-ph/990410( ) 
is generalized for any pairwise interaction. We describe two 
applications of the generalized framework. First, we obtain a 
class of protocols for selecting an arbitrary interaction term in 
an n-qubit Hamiltonia n. This class includes the scheme given 
in quant-ph/0106064v2. Second, we obtain a class of pro- 
tocols for inverting an arbitrary, possi bly unknown n-qu bit 
Hamiltonian, generalizing the result in quant-ph/0106085^1. 



I. INTRODUCTION 

An important element in quantum information pro- 
cessing is the ability to efficiently convert a set of primi- 
tives, determined by the physical system, to perform the 
desired task. In many physical systems, the primitives 
are "local manipulations" such as fast single qubit oper- 
ations that can easily be controlled, and a given nonlocal 
system Hamiltonian that cannot be changed. In this case, 
the desired task may be approximated or simulated by 
interspersing the given Hamiltonian evolution with local 
manipulations. The resources of simulation include the 
amount of local manipulations and the total operation 
time of the given Hamiltonian. 

Such simulation was extensively studied in the context 
of NMR quantum computation [|l]-[3| in which the natu- 
rally occurring Hamiltonian cannot be controlled. Refer- 
ence J2| presents a method based on Hadamard matrices 
to convert the evolution due to any Hamiltonian of the 
form 
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where a z 



to that of a particular term gi m <Jz 
is a Pauli matrix acting on the i-th qubit .n The same 
protocol applies universally for all coefficients gij and Wj. 
Other arbitrary evolutions can in turns be obtained by 



A similar method was reported independently in Ref. H. 



reduction to the universality construction j| || . The sim- 
ulation of the single term is exact, and does not require 
frequent local manipulations. A related task to stop the 
interaction is also addressed. The method aims at mini- 
mizing the amount of local operations, measured by the 
number of interval of free evolution or the number of sin- 
gle qubit gates required. 

A more general problem was addressed recently in 
Ref. M. A major step is to convert the Hamiltonian 
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the Pauli matrices. The single term is then used to sim- 
ulate the dynamics due to any other Hamiltonian 
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For each required accuracy level, both local and nonlocal 
resources are analyzed. 

Related problems were also discussed in Refs. ||||, fo- 
cusing on more specific types of given Hamiltonians and 
the operation time for the simulation. Bounds are de- 
rived for the operation time to use 
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to simulate an arbitrary Hamiltonian H' 2 given by Eq. (0). 
Methods and the required operating times for a Hamil- 
tonian H to simulate its inverse evolution e lHt are also 
studied]^] This is done for the Hamiltonians H% in Eq. (Q) 
and H4 given by 
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which is essentially the most general Hamiltonian for 
pairwise interaction with permutation symmetry.!] 



Throughout this paper, % — 1, and the time evolution due 



to a Hamiltonian H is given by e 



Note the — sign in 



the exponent. The inverse evolution is given by e l . This 
notation follows from the Schrodinger equation. 

3 The two-qubit case was independently considered in 
Ref. 



The general principle in these simulation schemes is 
to transform some coupling terms to the desired form 
and to cancel out the rest, by interspersing the free evo- 
lution with single qubit operations. In this paper, we 
generalize the framework developed in Ref. [Q to apply 
to the more general given Hamiltonian in Eq. (0). Using 
this framework, we find a class of schemes that select a 
term a\ ® o^ from a Hamiltonian given by Eq. (|2|). 
This class of schemes includes the method presented in 
Ref. jj. We also consider time reversal using the gen- 
eralized framework, and present a class of protocols to 
reverse an arbitrary Hamiltonian given by Eq. (0). They 
require an operation time c6nt of H2 to simulate the re- 
versed evolution e lH2t where 1 < c < 2 and c w 1 for 
large n. The schemes are designed to be independent of 
the given Hamiltonian, and are applicable even when H2 
is unknown. This generalizes the results in Ref. ||. 

After the initial submission of this paper, related work 
were independently reported in Ref. Jul. Methods to 
stop the evolution due to a Hamiltonian given by Eq. (||), 
and to select all diagonal coupling terms between a des- 
ignated pair of qubits are given. Simplifications for diag- 
onal couplings and higher order terms were considered. 
Our framework to stop the evolution is closely related 
to, but subtly different from that in Ref. |Ll|| , and allows 
more flexibility in selecting coupling terms, such as the 
selection of any single interaction term. We are also inter- 
ested in a tighter bound on the required local resources, 
and in the inversion of Hamiltonians. 

This paper is structured as follows. In Section y, we 
review the framework and various resulting schemes in 
Ref. 0], with a slight change from the original NMR 
based notations. The framework is generalized for any 
Hamiltonian given by Eq. (0) in Section [II. The first 
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application to select individual coupling terms from the 
given Hamiltonian is discussed in Section IV. The sec- 



ond application to simulate time reversal is discussed in 
Section M as a simple application. The technical details 
of the construction are described in Appendices |X| and 



II. SELECTIVE COUPLING USING HADAMARD 
MATRICES - A REVIEW 

A. Statement of the problem 

We review the method developed in Ref. . Consider 
an n-qubit system, evolved according to the Hamiltonian 



H Z =Y,9^ ] ®<^ 



(j) 

z > 



(6) 



i<j 



where gy are arbitrary coupling constants. The goal is 
to evolve the system according to only one term of the 
Hamiltonian: 



using single qubit operations. We call this task "selective 
coupling." This is closely related to the task of stopping 
the evolution or "decoupling." We first develop a frame- 
work for decoupling. Then we construct decoupling and 
selective coupling schemes using Hadamard matrices. 



B. Decoupling scheme for two qubits 

We motivate the general construction using the sim- 
plest example of decoupling two qubits. Let U t — e^ lH ' 1 , 
where the Hamiltonian is given by H z = 1712 to~ z (g) a z ■ 
We use the shorthand X^ 1 ' for ai ■ We also use the 
important identity 

Ue M tf = e UMuf (8) 

where M is any bounded square matrix and U is any uni- 
tary matrix of the same dimension. As the Pauli matrices 
anticommute, 

X {2) U t X {i) = X {2) e ~ i9l2t "P®"™ X^ 

_ e -igiat ^^(-of ) _ jj-i 

Thus adding the gate X^ before and after the free evo- 
lution reverses it, and 



X^U t X^ u t =i. 



(9) 



This illustrates how single qubit operations can trans- 
form the action of one Hamiltonian to another. 

Equation (|9|) can be written to highlight some essential 
features leading to decoupling: 

e -igilt(+a^)®{-af>) x e -ig 12 t{+a^)®(+a[ 2) ) _ MQ\ 

Each factor corresponds to a "time interval" of evolution. 

(i) 

1. In each interval, each o~ z acquires a — or + sign, 
according to whether X^ 1 ' are applied or not before and 
after the time interval. 

2. The bilinear coupling is unchanged (negated) when 
the signs of a z ' and a z ' agree (disagree) . 

3. Since the matrix exponents commute, negating the 
coupling for exactly half of the total time is necessary 
and sufficient to cancel out the coupling. 

The crucial point leading to decoupling is that, the 
signs of the a z matrices of the coupled qubits, controlled 
by the X gates, disagree for half of the total time elapsed. 



C. Sign matrix and decoupling criteria 

We now generalize the framework for decoupling to n 
qubits. Each of our schemes concatenates some m equal- 
time intervals. In each time interval, the sign of each 



(Tz can be + or — as controlled by the X gates. Each 
scheme is specified by an n x m "sign matrix" , with the 
(?, a) entry being the sign of <7z in the a-th time inter- 
val. The entries in each column represent the signs of 
all the qubits at the corresponding time interval and the 
entries in each row represents the time sequence of signs 
for the corresponding qubit. We denote a sign matrix for 
n qubits by S n - For example, the scheme in Eq. (110) can 
be represented by the sign matrix 



So 
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Following the discussion in Sec. [IB, we have 



Decoupling criteria version 1 Decou- 
pling is achieved if any two rows in the sign 
matrix disagree in exactly half of the entries. 

For example, the following sign matrix decouples four 
qubits: 



Si = 



+ 
+ 



(12) 



More explicit, the scheme is given byn 

U t x (X^X^ UtX^X^) x (13) 

(X^X^U t X^X^) x (X^X^U t X^X^), 



where Ut 



-iH z t 



has six possible coupling terms in the 



4-qubit Hamiltonian. The relation between the scheme 
and S4 is illustrated in Figure pi. 
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FIG. 1. (a) Converting the sign matrix S* to the scheme 
in Eq. (|13[). A "— " sign in the i-th row and a-th column 
translates to X^'' (acting on the i-th qubit) before and af- 
ter the a-th time interval, (b) Simplifying the scheme using 



4 Note that the commuting factors in Eq. (h3|) are arranged 
to visually correspond to the sign matrix. 



From now on, we only consider the sign matrices, which 
completely represent the corresponding schemes. For n 
qubits when n is large, n X m sign matrices with small 
m can be difficult to construct. For example, 



S n 



+ 



+ 



in which an interval from a previous row is bifurcated 
takes m = 2 n_1 . As the number of columns represents 
local resources for the simulation, our goal is to find so- 
lutions with the smallest number of columns. We solve 
this problem by first rephrasing the decoupling criteria: 

Decoupling criteria version 2 Identify- 
ing ± with ±1 in S n , decoupling is achieved 
if any two rows have zero inner product, or 

We now present very efficient solutions to the decoupling 
criteria, which are derived from the Hadamard matri- 
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D. Hadamard matrices and decoupling scheme 

A Hadamard matrix of order n, denoted by H(n), is 
an n x n matrix with entries ±1, such that 



H(n)H(n) T = nl . 



(14) 



Thus every H{n) 1 if exists, is a valid sign matrix corre- 
sponding to a decoupling scheme for n qubits using only 
n time intervals. The following is a list of interesting facts 
about Hadamard matrices (see Refs. |12|-fl5| for details 
and proofs). 

1. Equivalence Any permutation, or negation of any 
row or column of a Hadamard matrix preserves the 
orthogonality condition. Thus each Hadamard ma- 
trix can be transformed to a normalized one, which 
has only +'s in the first row and column. 

2. Necessary conditions H(n) exists only for n = 1, 
n = 2 or n = mod 4. 

3. Hadamard's conjecture u& H{n) exists for every 
n = mod 4. This famous conjecture is verified for 
all n < 428. 
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If H(n) and H{m) 
H(m) is a possible H(nm). In 



4. Sylvester's construction 
exist, then H(n) 
particular, H(2 r ) can be constructed as H(2)® r , 
which is proportional to the matrix representation 
of the Hadamard transformation for r qubits. 



5. Paley's construction [|18| Let q be an odd prime 
power. If q = 3 mod 4, then H(q+l) exists; if 
q = 1 mod 4, then H(2(q+1)) exists. 

6. Numerical facts p2[ For an arbitrary integer n, 
let n be the smallest integer satisfying n <n with 
known H(n). For n < 1000, H(n) is known for 
every but 6 possible orders, and n — n < 7. For 
n < 10000, if (n) is unknown for only 192 possible 
orders and n — n < 31. 

The nontrivial existence of so many Hadamard matrices 
may be better appreciated by examining the following 
example of H (12), obtained with Paley's construction. 
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Thus there is a simple decoupling scheme for n qubits if 
H(n) exists. Using Hadamard matrices, decoupling and 
recoupling schemes for an arbitrary number of qubits can 
be easily constructed, as will be shown next. 



E. Decoupling and selective coupling 

Decoupling When an H(n) exists, it corresponds 
to a decoupling scheme for n qubits concatenating only 
n time intervals. When an H(n) does not necessarily 
exist, consider H (n) and choose any n rows to form an 
S n - Then S n corresponds to a decoupling scheme for n 
qubits requiring n time intervals. As an example, Sg can 
be chosen to be the first nine rows of H(12) in Eq. (15): 
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Selective coupling To implement selective coupling 
between the i-ih and the j-th qubits, any two rows in the 
sign matrix should be orthogonal, except for the i-th and 



j-th rows that are identical. The coupling 
acts all the time while all other couplings are canceled. 
The sign matrix can be obtained by taking n — 1 rows 
from H(n—1). For example, to couple the last two among 
9 qubits, we can take S% to be the 9x8 matrix obtained 
from appending the last row of H(8) to itself. Alterna- 
tively we can take the 2-nd to the 9-th rows of -ff (12) in 
Eq. (jig), and repeat the last row: 
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The extra feature of this Sg is that, all row sums are zero. 
This is because -ff(12) in Eq. ( ^5[ ) is normalized, so that 
all rows except for the first have zero row sums. This 
automatically removes any local (linear) terms J ^2 i LO{Oz 
in the Hamiltonian, without extra local resources (see 
Ref. Q for a full discussion). Finally, we note that cou- 
pling terms involving disjoint pairs of qubits can be se- 
lected simultaneously. 



F. Discussion 

Upper bound on n For n qubits, selective coupling 
requires at most n intervals and nn single qubit gates. In 
fact, n = en where c is very close to the ideal lower bound 
c = 1. First, if Hadamard's conjecture is proven, n only 
depends on n mod 4, and V n n—n < 3. Even without this 
conjecture, the present knowledge in Hadamard matrices 
implies n-n < 8 V n < 1000, and n-n < 32 V n < 10000. 
A detailed proof for c w 1 for arbitrarily large n is given 
in an Appendix of Ref. pi , while Sylvester's construction 
puts an immediate loose bound of c < 2. 

Gate simulation vs dynamics simulation pi The 

previous discussion assumes that the goal is to simulate 
the final unitary transformation due to the Hamiltonian 

(i) (i) 

o~z <8> <Jz for time t. Due to the commutivity of all the 
possible coupling terms, we only need to divide the time 
into n time intervals, each with finite duration t/n. On 
the other hand, if the goal is to simulate the dynamics 
due to a z (Ssa z for time t, one should instead apply the 
scheme to simulate the unitary gate "e - " 7 ** ® CT ^ At " t/At 
times where At is a small time interval. We will focus on 
simulating the dynamics of a Hamiltonian in the rest of 
the paper. 



III. GENERALIZED FRAMEWORK FOR 
ARBITRARY iV-QUBIT HAMILTONIANS 

We now generalize the framework for a more general 
riven Hamiltonian (Eq. (0)): 



product S x .*S y = S z . The (i, a) entry of S a 
is the sign of o~a in the a-th time interval. 



H2 = Y J Ha^p^ ) 
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The goal is again to simulate the evolution due to one 
specific coupling term o~\ ' (g) oVj . Passing from Eq. (|l|) 
to Eq. (|l6|), the first difference is the noncommutivity 
of the terms in Eq. (|lq). The second difference is the 
presence of all three Pauli matrices acting on the same 
qubit, besides a much larger number of coupling terms. 
We adopt a common approach [|7|,p|JlC|l that employs 
sufficiently frequent local manipulations to make the ef- 
fect of the non-commutivity negligible. This is based on 
the identity 



We omit the number of qubits, n, in S x , S y , S z for 
simplicity of notation. The entry-wise product .* of 
two matrices is also known as the Schur product or the 
Hadamard product. 



IV. SELECTIVE COUPLING FOR N QUBITS 
WITH ARBITRARY PAIRWISE COUPLING 

Under the generalized framework, we state the criteria 
for decoupling and selective coupling for n qubits: 

Criteria for decoupling and selective 
coupling Decoupling is achieved if any two 
rows taken from S x , S y , S z are orthogonal. 
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which implies that effects of non-commutivity are of 
higher order in the small U in Eq. ([17]) . Thus the discus- 
sion proceeds neglecting the non-commutivity, and prod- 
ucts of unitary evolutions are replaced by sums of the 
exponents. With this simplification, the framework in 
the previous section is readily generalized. 

Again, we consider a class of schemes that concatenate 
(short) equal time intervals of evolution. The essential 
features for decoupling are as before: 
1. In each interval, each a a acquires a + or — sign, 
which is controlled by the applied local unitaries to be 
described. 
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Selective coupling of ct 7 <g> aij is achieved if 
the l-th. row of 5 7 is identical to the fc-th row 
of S v , but any other pair of rows from S x , S y , 
S z are orthogonal. Local terms are removed 
if all row sums are zero. 

Sign matrices S XtVjZ satisfying the criteria can be con- 
structed from special Hadamard matrices endowed with 
certain extra structures. We now describe the construc- 
tions of the sign matrices, which elicit the special prop- 
erties required of the starting Hadamard matrices. The 
more difficult and technical constructions of these special 
Hadamard matrices are given in detail in Appendices |A| 
and |B|. 

Suppose we want to decouple n qubits using a 
Hadamard matrix H(m). The orthogonality condition is 
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changed (negated) when the signs of a a and a^ ' agree 
(disagree) . 

3. To the lowest order in the duration of the time inter- 
vals, negating a coupling for exactly half of the intervals 
cancels it. 



2. The bilinear coupling gij a p era ® erg for i ^ j is un- automatically satisfied if the rows of S x , y ^ z are taken to be 

distinct rows of H(m). It remains to ensure S x .*S y = S z . 
We call a set of 3 vectors of equal length and with entries 
±1 a Schur-set if they entry- wise multiply to + + • • • +. 

For example, {[H ], [ — I — ], [ h]} is a Schur-set. 

If H(m) has at least 3n rows that partition into n Schur- 
subsets, one can ensure S x .*S y = S z by choosing the z-th 
rows of S x , S y , S z to be the rows of the i-th Schur-subset. 
This poses the first extra property on H(m) - its rows 
partition into many Schur-subsets. Note the immediate 
lower bound on the size of H(m), \_(m — 1)/3J > n under 
this construction. This is not a strict lower bound, as a 
modification (to be described later) to the construction 
enables n to be replaced by n + 1 , but the modification 
removes some useful properties. 



The 
one 



a=x,y,z 



generalized framework differs from the original 
I in that, the signs of the three Pauli matrices 
acting on the same qubit i in the same time 



interval are not independent. In fact, the three signs mul- 
tiply to +, because conjugating ?(*' ■ <?M by (local) uni- 
taries on the z-qubit transforms r "' by an SO (3) matrix. 
Conversely, any sign assignment satisfying this constraint 
can be realized. The possible signs for a x ,a y ,a z are 



+ + +, H , — I — , h, and are realized by applying 

/W, a x , a y , o~z respectively before and after the in- 
terval. Incorporating these considerations, we generalize 
the previous framework: 

A scheme for n qubits that concatenates m 
intervals can be specified by three n x m sign 
matrices S x , S y , S z , related by the entry-wise 



Schemes for selective coupling can be derived from de- 
coupling schemes as follows. Let {7,77, v\ = {x, y,z} be 
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distinct labels. To select the coupling <r 7 
modify the decoupling sign matrices as follows: 
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l-th row 
k-th row 



l-th row 
k-th row 

where v a 



Ury .* Uni 



for a 6 {7, 77, v} denote the rows. To select 
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the coupling a\ ' ® <r 7 , proceed as before and further 
swap the fc-th rows of 5 7 and S v (these are w 1 and v 1 
in the lower diagram). In both cases, it is necessary to 
ensure w 7 .*w 7 is a row of H{rn) and is not used elsewhere 
in S x , y ,z- To ensure w 1 .*v 1 is a row of H(m), one poses a 
second special property in H(m) that, it contains 5 rows 
/1 ... 5 such that 



/l-*/2 — f3-*fi — fl 



5 • 



(18) 



Then we can simply choose u 7 = /s, u^ = /i ( 2 in the 
decoupling scheme, and replace u> 7 , u 7 .*u; 7 by fe^. To 
ensure /^ are not used elsewhere in S x ,y,zi the simplest 
method is to exclude the Schur-subsets originally contain- 
ing them. This is often unnecessary, as /1 ... 5 can often 
be found such that fy^ are not in any Schur-subset in the 
starting decoupling scheme. Meanwhile, the w XjV ^ z can 
be "recycled" for an extra qubit. Altogether, the scheme 
can handle nil qubits depending on the situation. 

The above constructions of decoupling and selective 
coupling schemes involve only rows with zero sums if the 
Hadamard matrix H(m) is normalized. This automati- 
cally removes all local (linear) terms without extra local 
manipulations. It also means that, one can append an 
extra row + + •••+ to all of S XtVtZ to handle an extra 
qubit, the local terms of which, if exist, have to be dealt 
with outside of the scheme. 

We now summary the results in Appendices and 
H concerning the difficult issue of constructing normal- 
ized Hadamard matrices that consist of many Schur- 
subsets and have 5 rows satisfying Eq. (jig)- The sim- 
plest of these are the Sylvester-type Hadamard matri- 
ces H(2)® r (referred to as Sylvester matrices from now 
on). In Appendix EL we show how H{2)® r can be par- 
titioned into (2 r — l)/3 and (2 r — 5)/3 Schur-subsets 
when r is even and odd respectively. When construct- 
ing a scheme for n qubits based on Sylvester matrices, 
r is chosen to satisfy the inequalities which are approxi- 



mately^] 2 r_1 < 3n < 2 r , and the number of intervals is 
approximately between 3n — 6n. Appendix M describes a 
more involved construction that combines a normalized 
generalized Hadamard matrix with a Hadamard matrix 
having all the required properties to produce a larger 
Hadamard matrix having the same properties. This sig- 
nificantly improves on the worst case number of intervals. 
For asymptotically large n, the number of intervals, 3cn, 
has c w 1. The value of c, as a function of n, is plotted 
in Figure @. 
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FIG. 2. A plot of c as a function of n, where 3cn is the 
number of intervals for an n-qubit scheme. The lower bound 
of c is 1. Using the Sylvester construction alone, c fluctuates 
between 1 and 2 indefinitely. 

In many applications, it may be useful to select more 
than one coupling term from the Hamiltonian. For ex- 
ample, one may select all coupling terms between the i-th 
and the j'-th qubits in H2, by choosing the i-th and the 
j'-th rows of S x ^ v>z to be identical, and any other pairs of 
rows to be orthogonal. These 6 identical rows have to be 
+ + •••+ in order to satisfy S x .*S y — S z . The required 
number of intervals is the same as decoupling n— 2 qubits. 
This pairwise coupling can further be changed to a de- 
sired form following methods in Refs. Jl^Q. Depending 
on circumstances, this composite method can reduce the 
operation time of the given Hamiltonian by a constant 
factor but can potentially increase the local resources by 
a constant factor. As a final remark, the selective cou- 
pling schemes in both Refs. 0jll| are to first select the 
diagonal coupling terms or all coupling terms between 



5 The exact number depends on whether one is concerned 
with decoupling, selective coupling, or inversion of Hamiltoni- 
ans, and whether local terms are to be handled by the scheme. 
These may affect n by a difference of or ±1. To avoid the 
cumbersome description of all possible variations, from now 
on, we simply give approximates and omit ±1. The readers 
can work out the exact bounds tailored for their tasks. 



a pair of qubits before further manipulations. The ini- 
tial selection is done essentially using some version of 
the Sylvester construction, with even and arbitrary r in 
Refs. and |ll[ respectively. 

V. UNIVERSAL TIME REVERSAL 

We now apply the framework of sign matrices to simu- 
late time reversal. As a first example, we apply the orig- 
inal framework in Ref. B to reverse H i given by Eq. (TO) . 
First construct a decoupling sign matrix S n for n qubits 
using a normalized Hadamard matrix H(n+1), excluding 
its first row. In this case, all entries in the first column 
of S n are "+" . A sign matrix for time reversal, denoted 
by T n , is obtained by removing the first column of S n . 
In S n , any two rows have zero inner product, thus any 
two rows in T n have inner product —1, and any coupling 
term is reversed by exactly the same amount (note the 
significance of the — sign in the inner product). Further- 
more, each row of S n has zero row sum, so, each row 
in T n has row sum —1, and the local terms are reversed 
by the same amount as well. Therefore, any Hamiltonian 
given by Eq. (FA) can be reversed. The simulation requires 
ni = Ti+1 — 1 intervals, and simulates the reversal for the 
duration of one interval. Hence the simulation factor |1C| ] 
is 1/nj or the overhead Ppl is nj. The overhead is be- 
tween n and n + 3 if Hadamard's conjecture is true, and is 
en for c ~ 1 in any case. Without local terms, nj = n— 1 
ranges between n — 1 and n + 2 if Hadamard's conjecture 
is true. Note that the protocol is independent of Hi, and 
can be applied to any Hi which can even be unknown. 
We remark that a lower bound n — 1 was derived for 
the much more specific and known Hamiltonian H$ in 
Eq- @j |J- In view of this, the generality of the present 
result comes almost for "free" . 

To reverse a Hamiltonian given by Eq. (J2J) using the 
generalized framework, we again construct sign matrices 
Sx,y,z for a decoupling scheme (Section [TV]) with zero 
row sums. Furthermore, our constructions can be made 
to ensure all entries in the first columns of S x ,y,z are +• 
Therefore, removing the first columns in S x>VtZ results 
in a time reversal scheme. The overhead is again the 
number of intervals in the reversal scheme, which is 3cn 
for c rs 1 for large n. Again, the same protocol applies 
to any Hi and thus applies to unknown Hamiltonians as 
well. 

As a comparison, the reversal method reported in 
Ref. M for H4 given by Eq. (H) depends on the knowledge 
of d a , and whether they have the same sign or not. When 
all d a have the same sign, the overhead is maxirf 1 — 
— - — u — — which ranges from w n — in. The gen- 



or 



max \d a 



knowledge of the given Hamiltonian. 



VI. CONCLUSION 

We have generalized the framework for Hamiltonian 
simulation and the methods for decoupling and selec- 
tive coupling in Ref. 0. We rederive, as a special case 
of our construction, the crucial step of selecting a cou- 
pling term in the simulation of n-qubit Hamiltonians in 
Ref. [Q . We also apply the technique to extend the time 
reversal problem considered in Ref. J9| from permutation 
invariant purely nonlocal Hamiltonians to an arbitrary 
n-qubit Hamiltonian. 

Our framework based on sign matrices allows the 
complicated criteria for various simulation tasks to be 
rephrased in relatively simple orthogonality conditions, 
for which solutions can be obtained with the connections 
to Hadamard matrices. 
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APPENDIX A: SCHUR-SUBSETS IN 
SYLVESTER MATRICES 

In this Appendix, we study the properties of the 
Sylvester matrices useful for decoupling, selective cou- 
pling or inversion of Hamiltonian. First of all, they are al- 
ways normalized. We now describe a method to partition 
the rows of the Sylvester matrix H(2)® r into (2 r - l)/3 
and (2 r — 5)/3 Schur-subsets when r is even and odd re- 
spectively. It will be obvious from the construction that 
each H(2)® r contains 5 rows satisfying Eq. (jlq). 

First, we introduce some notations. Let {0, 1} be the 

index set for the rows and columns of H (2) = 

For example, the (0, 0) entry is 1. We use the shorthand 
Hij for the (i, j) entry of H(2). Therefore Hy = (-l) ij '. 
Likewise one can label the rows and columns of H(2)® r 
with composite indices i = («'i,i2, • ■ ■ , v) which are r-bit 
strings. We have 

"ij = "iiji x -"i2j2 X ' ' ' X "irjr 

= f ]_"\*lJl+*3jaH Mrjr 

= (-l) ij (Al) 

where i • j denotes the usual inner product of i and j. For 
each 1, 

H n x Hji x Hid = (_i)«+J-i+k-i = (_!)(i+J+k).i ( A2 ) 



Therefore, the i-th, j-th, k-th rows form a Schur-set iff 
Vi H n x Hji x iJ kl = 1 iff i © j © k = 0. We refer to such 
a triple of r-bit strings as a Schur-set also. The problem 
reduces to showing that, the set of non-zero r-bit strings 
partitions into (2 r — l)/3 and (2 r — 5)/3 Schur-subsets 
when r is even and odd respectively. This can be proved 
by separate inductions on the even and odd values of r. 

For even values of r > 2, the induction hypothesis 
(IH) is that, the set of all r-bit strings partitions into 
K = (2 r - l)/3 Schur-subsets Eh, E 2 , ■ ■ •, E K , and the 
singleton {0}. The IH is clearly true when r = 2. Sup- 
pose it is true for some even r > 2. For each Schur-set 
Ei = {*i, *2j ^3} of r-bit strings, we can obtain 4 Schur- 
sets of (r+2)-bit strings: 

{Olti, 10*2,11*3} 
{0lt 2) l0i3 S llii} 
{Oils, lOii.llta } 
{00«'i,00i 2 ,00i 3 } (A3) 

We also have an addition Schur-set {010, 100, 110}. Al- 
together, we find 4(2 r - l)/3 + 1 = (2 r + 2 - l)/3 Schur- 
sets of (r+2)-bit strings that include all strings except for 
000. This completes the induction when r is even. For 
odd values of r > 3, the IH is: 

1. the set of all r-bit strings partitions into K = (2 r — 
5)/3 Schur-subsets E\, E 2 , ■ ■ ■, Ek, and a set of 5 
remainders F = {wi, w 2 , W3, W4, 0}, 

2. Ek — {ki,k2,k 3 } satisfies 

k\ = W\ © W2 = w 3 © W4 , 

k 2 = W\ © w 3 = w 2 © wa . (A4) 

For r = 3, the IH is true, for example, by putting E\—k = 
{001, 100, 101}, and F = {010, Oil, 110, 111,0}. Suppose 
the IH is true for some odd r > 3. Using the method of 
Eq. (O), we obtain A(K - 1) Schur-subsets of (r + 2)- 
bit strings from Eh., ■ • •, Ek-i- There are 32 remaining 
strings, {00,01, 10, 11} x {fe 1: fc 2 , k 3 ,wi,W2,W3,W4,0}. If 
we represent distinct bit strings as distinct points, and 
Schur-sets as triangles, the IH implies the following rela- 
tions for the 6 r-bit strings: 





(A5) 



Note that two triangles that share a common vertice 
represent two Schur-sets that are not disjoint. From 



Eq. (A5), we can easily find another 8 disjoint Schur- 
subsets of (r+2)-bit strings: 



lluii lltua 




01t«i 00u; 2 




llu>3 QlWi OOu/3 11^4 

OOtoi 01u) 2 

IOwji 10w 2 




01l(J 3 



10w 4 10«j 3 




00104 



(A6) 



Denote the remaining 8 strings as k[ — 110, k' 2 — Oiks, 
10fc 3 , tt>i = 00fc 3 , w' = llfc 3 , w' 3 = 01O, w' = 100, 



A'., 



0' = 00O. An additional Schur-set {/e^fc^/cg} can be 
formed, and the 5 remaining rows w' i= 



,1,2,3,4' 



0' 



not form any more Schur-set. Moreover, k'^w^ satisfy 
Eq. (A4). Finally, the number of (r+2)-bit Schur-sets is 

(2 r + 2 



K' = 4((2 r - 5)/3 
the induction. 
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5)/3, completing 



The inductive proof, together with the identification 
Hy = (— 1)*'J, provides a constructive method to parti- 
tion the rows of H (2)® r into Schur-subsets. 

Finally, the fact H® r contains 5 rows satisfying 
Eq. (|J) is manifest for odd r in Eq. (|Aj). In fact Eq. (|l|) 
is exactly equivalent to a pair of connected triangles as 
in Eq. flA5| ). The case for even r follows by adding a "0" 
to each point in the joint triangles for odd r in Eq. (|A5|). 



APPENDIX B: CONSTRUCTION USING 
GENERALIZED HADAMARD MATRICES 

In this Appendix, we construct Hadamard matrices of 
orders other than 2 r with the properties required of our 
schemes. This is done by composing a Hadamard matrix 
with a generalized Hadamard matrix |12f| . 

Let G be a group of order g, with group operation *. A 
generalized Hadamard matrix over G 12,13], GH(g,A), 
is a gX x gX matrix whose entries 7^ are elements of G, 
and for i ^ j, the sequence {7^ * 7jj : 1 < I < .9A} 
contains each element of G X times. This sequence is the 
entry-wise "division" of the i-th row by the j-th row. For 
example, a Hadamard matrix H(4n) is a GH(2, 2n) over 
the multiplicative group {+1,-1}. 

We are interested in generalized Hadamard matrices 
over GF(4), with g = 4 elements written in an unusual 
manner: 



+ + 

+ ,-,+ ,- 

+ - + 



(Bl) 



The group operation .* is the entry-wise multiplication 
for the triples. The GH(4, A) is a 4A x 4A array of triples 
(written as a column vector). For example, a possible 
GH(4, 1) is given by 



1 = 





+ 
+ 
+ 

+ 
+ 
+ 

+ 
+ 
+ 

+ 
+ 
+ 




+ 
+ 
+ 

+ 
+ 

+ 




+ 
+ 
+ 

+ 

+ 
+ 




+ 
+ 

+ 

+ 
+ 

+ 





(B2) 



We state some useful facts about generalized Hadamard 
matrices (see for example, Ref. [[L2[ Part IV). 

1. Equivalence pl| Permuting rows or columns, or 
multiplying any row or column by a fixed element 
from the center of G preserve the defining proper- 
ties of a GH(g, A). Therefore, any GH(g, A) over an 
abelian group G is equivalent to a normalized one, 
in which all entries in the first row and column are 
equal to the identity element in G. 

2. Kronecker product |20| , |21| The Kronecker product^] 
of GH(g, Ai) and GH(g, A2) over the same abelian 
group G is a GH(g, 5A1A2) over G. 

3. [£2| Let q = gX — 1 be a prime power. If GH(g, A) 
over G exists, Vj>o GH(g, Xq l ) over G exists. 

Applying the above facts to GF(4), the existence of 
GH(4,A) for A = 1,2, implies that for A = 1, 2, 3, 4, 
8, 9, 12, 14, 16, 24, 27, 33, ■•■. 

Let H(m) be a Hadamard matrix with all the re- 
quired properties: normalized, having 5 rows fi satisfy- 
ing Eq. (18), and having some 3n rows forming n disjoint 
Schur-sets. We can assume that the rows of H(m) are 
ordered so that members of each Schur-subset occur con- 
secutively, starting from the first row, and the last m — 3n 
rows do not belong to any Schur-subsets. For example, 
the 1-st to 3-rd rows form a Schur-subset, and same for 
the 4-th to 6-th rows, and so on. Converting every 3 con- 
secutive rows into a row of triples. Do this for the first 
3n rows of H(m). The resulting array He is an n x m 
array of triples. Call the last m — 3n rows of H(m) Hp- 
For example, 



6 The Kronecker product in this case is defined with the 
group operation * replacing the usual multiplication. 



H(2) 



■Si 2 



H h 

+ H 

H 



H E = 

H F =[+ + + +} 



(B3) 



As a second example, consider H® A . The correspond- 
ing Hp is [+ • • • +]. Following the previous Appendix, 
Schur-subsets of 4-bit strings are {0101,1010,1111} 
{0110,1011,1101} {0111,1001,1110} {0001,0010,0011} 

-(0100, 1000, 1100}. Thus, H E corresponding to H(2) m 
is equal to 



+ 


— 


+ 


— 


+ 


— 


+ 


— 


+ 


— 


+ 


— 


+ 


_ 


+ 


— 


+ 


+ 






+ 


+ 






+ 


+ 






+ 


+ 






+ 


— 


— 


+ 


+ 


— 


— 


+ 


+ 


— 


— 


+ 


+ 


— 


— 


+ 


+ 


— 


+ 


— 


— 


+ 


— 


+ 


+ 


— 


+ 


— 


— 


+ 


_ 


+ 


+ 


+ 


- 
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+ 


- 


- 


- 


- 
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+ 
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— 
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+ 
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+ 
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+ 


















+ 


+ 


+ 


+ 



Returning to the general construction, we now com- 
pose a normalized 7 = GH(<7, A) and H(m) to form a 
new Hadamard matrix: 

1. Take the Kronecker product He ® 7 (under .*). 

2. Convert each row of triples in He <8> 7 back into 3 
rows of ± (an operation we call "leveling"). Call 
the resulting matrix H' E . It is 3ra4A x m4A. 

3. Take the first coordinates of each triple in 7 and 
form a 4A x 4A matrix of entries ±. Take the (usual) 
Kronecker product of Hp with the above matrix to 
obtain H' F . It is (m — 3n)4A x m4A. 

4. Append H' F to H' E to obtain an m4A x ra4A matrix 
of ±1. Call this H'. 

We first show that the resulting matrix is indeed a 
Hadamard matrix. Denote by e^ the (i,j) entry of the 
leveled 7, and hij the (i,j) entry of H(jn). Then, H' E is 
explicitly given by: 





hu 

tl21 
Jl31_ 


.* 


7 






h%2 

fl22 
Jl32_ 


.* 


7 








ha 

h 5 i 

fl61 


.* 


7 






h-42 
h 5 2 
h%2 


.* 


7 





























(B4) 



or 



" hu £11 hu £12 ■ ■ • 

^21 £21 hll £22 ' ■ ■ 
hzi £31 hzi £32 

hll £41 hn £42 ■ ■ ■ 
ft-21 £51 ^21 £52 ' ■ ■ 
^31 £61 hu £62 ' ■ ■ 


^12 £11 ?»12 £12 ' • ' 
^22 £21 h22 £22 • • • 
^32 ^31 ^32 £32 ■ ■ ' 

^12 £41 hl2 £42 • • • 
^22 £51 ^22 £52 ' ■ ' 
^32 £61 ^32 £62 • • • 




/l41 £n /141 £12 • ■ ■ 
/l51 £21 /l51 £22 ' ■ ■ 
^61 £31 h(,i £32 ■ ■ ■ 

h.41 £41 hu £42 ■ ■ ■ 

hisi £51 hsi £52 • ■ ■ 

/l61 £61 ^61 £62 ■ ■ ■ 


^42 £11 ?»42 £12 ' ■ ' 
^52 £21 ^52 £22 • • • 
^62 £31 ^62 £32 ■ ■ ■ 

h.4,2 £41 ft-42 £42 • • • 
h$2 £51 h<S2 £52 ■ ■ • 
^62 £61 ^62 £62 • • • 











(B5) 



From Eq. 



the rows of H' E are of the form: 



[hn,h i2 ,- •■ ] <8> [eji,ej2, 



(B6) 



where i e { 1 , • 



The rows of H' F 



in}, j e {1, • • • , 12A} and i — j mod 3. 

are of the same form, with i G {3n + 
1, • • • , m}, j = 3fc + 1 for k £ {0, • • • , 4A - 1}. Thus 
each row in H' is specified by i and j. Note that if i j^ i\ 
[ha,hi2, • • • ] and [hi>i,hi>2, ■■■ ] are orthogonal. If j ^ j' 
and j = f mod 3, then [eji, 6j2, • • • ] and [eyi, Cj'2, • • • ] 
are orthogonal since 7 is a generalized Hadamard matrix. 
Therefore, any two distinct rows in H' F are orthogonal 
(orthogonal in at least one tensor component), and any 
row from H' E is orthogonal to any row from H' F (orthog- 
onal in the first tensor component). Finally, let i,j and 
i',j' specify two distinct rows in H' E . If i ^ i', they are 
orthogonal in the first tensor component. If i = i 1 , then 
j 7^ fi j = j' mod 3, and they are orthogonal in the 
second tensor component. Thus H' has orthogonal rows 
and is a Hadamard matrix. f\ 

Note that the rows of H' E completely partition into 



7 Readers who are familiar with combinatorial designs may 
interpret the core part of the above construction as taking the 
Kronecker product of an orthogonal array (He) and a differ- 
ence matrix ( 7) , both over GF(4), to obtain a larger orthogo- 
nal array [ |12|Jl9| ] . The above construction is more elementary 
and allows slightly more features. 
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Schur-subsets due to the group structure of the triples 
under .*. The ratio of the maximum number of qubits 
handled by the scheme to the number of intervals is the 
same in H' and H(m) - everything is rescaled by a factor 
of 4A. This construction is therefore most useful when 
H(m) has a large fraction of rows forming Schur-subsets, 
and when A is not a multiple of 4. 

One can verify that the above construction results in a 
normalized H' up to a permutation of the rows. Further- 
more, because 7 is normalized [ha, hi2, ■ ■ ■ ]<&[+ + ■■■ ] 
occurs in H' for all i, Therefore, H' contains 5 rows 
fi ® [+ H ] that satisfy Eq. @. 

Let the number of intervals required for an n-qubit 
scheme be 3cn, and consider c as a function of n. We now 
use the above construction to put a loose upper bound 
on c for large n. A tighter bound on c for smaller values 
of n is plotted in Figure g. A scheme formed by com- 
posing H(m) = H{2)^ r -^ and 7 = GH(4,3<* +1 )) has 
2 r 3 (t+1) intervals, and handles at least (2 r -20)3* qubits. 
For a given n, we find the smallest value of (2 r — 20)3* 
larger than n. Let n = 2 r ° +A where A £ [0,1), and 
e = log 2 (27(2 r - 20)). We look for the smallest non- 
negative value of 

log 2 ((2 r -20)37n) = r-e + ilog 2 3-r - A (B7) 

where t, r are positive integers. When n is large, so is 
r . Since log 2 3 is irrational, {(tlog 2 3) mod 1} is a dense 
subset of [0, 1). We can find some t < r /21og 2 3 such 
that ilog 2 3 = N + S, S G [0, 1), and 6 - (A + e mod 1) 
is small and nonnegative. We also choose r — r — N + 
[A + e\ . Note that r is at least r /2. Then, 

log 2 c = log 2 (2'-3 t + 1 /3n) 

= r + t log 2 3 - r - A 

= r -A^+LA + eJ+^ + (5-r -A 

= [A + e\ + S - A - e + e 

= [8 - (A + e mod 1)] + e (B8) 

Both terms in the last line can be made small when r is 
large, and c«l. 
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